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ABSTRACT 

One mechanism suggested for the L to T dwarf spectral type transition is 
the appearance of relatively cloud-free regions across the disk of brown dwarfs 
as they cool. The existence of partly cloudy regions has been supported by 
evidence for variability in dwarfs in the late L to early T spectral range, but no 
self-consistent atmosphere models of such partly cloudy objects have yet been 
constructed. Here we present a new approach for consistently modeling partly 
cloudy brown dwarfs and giant planets. We find that even a small fraction 
of cloud holes dramatically alter the atmospheric thermal profile, spectra, and 
photometric colors of a given object. With decreasing cloudiness objects briskly 
become bluer in J — K and brighten in J band, as is observed at the L/T 
transition. Model spectra of partly cloudy objects are similar to our models 
with globally homogenous, but thinner, clouds. Hence spectra alone may not be 
sufficient to distinguish partial cloudiness although variability and polarization 
measurements are potential observational signatures. Finally we note that partial 
cloud cover may be an alternative explanation for the blue L dwarfs. 

Subject headings: brown dwarfs — stars: atmospheres 



-2- 



INTRODUCTION 



As brown dwarfs cool over time, their atmospheres undergo a sequence of chemical and 
physical changes that result in an evolving emergent spectrum and-consequently-varying 
spectral types. The most remarkable changes take place at the L to T type transition. 
The latest L dwarfs have red near-infrared colors, strong CO absorption in K band and 
relatively shallow water absorption bands modulating their spectra. Over a small effective 
temperature range of only 100 to 200 K the spectrum rapidly changes to exhibit blue near- 
infrared co lor, weakening CO absorption, strengthening CH4 absorption, and deeper water 
bands (see iKirkpatrick I (120051 ) for a review) . 



The proximal cause of these changes is the loss of cloud opacity. As clouds dissipate 
the visible atmosphere cools, bringing on the chemical change from CH4 to CO. Without 
clouds providing a significant, nearly gray opacity, flux can emerge through molecular opacity 
windows in the J and H bands with a brightening leading to a blueward color shift (Dahn 
et al. 2010, Tinney et al. 2003, Vrba et al. 2004). Studies of L + T binary dwarfs (e.g., 
Liu et al. 2006; Looper et al. 2008; Stumpf et al. 2010 and references therein) show that 
this brightening is an intrinsic signature o f the transition and not the manifestation of some 
other effect in color- magnitude diagrams ( Stumpf et al. 201ol ). 



Two main underlying causes of this loss in cloud opacity have been suggested. In 
one view the atmospheric dynamical state changes, resulting in larger particle sizes that 
more r apidly 'rain out' of the atmosphere, leading to a sudden clear ing or collapse of the 

This view is sup- 



cloud ( IKnapp et al.l l2004t iTsuji &: Nakajima 



2003 



Tsuii et al. 



2004 



ported by fits of spe c tra to model spectra (jSaumon fc Marley 1 120081 ) computed with the 



Ackerman fc Marley I ( 1200 ll ) cloud model. In that formalism a tunable parameter, f se< i con- 
trols cloud particle sizes and optical de pth. Larger / Sf d yield s lar ger particles along with 



physically and optically thinner clouds. ICushing et al.l ( 120081 ) and IStephens et al. I (120091 ) 



have demonstrated that progressively later T dwarfs can be fit by increasing / se( j across 
the transition. Late T dwarfs are generally best fit by models that neglect cloud opacity. 
Burrows et al.l ( 120061 ) also suggest changes in cloud particle sizes as a possible mechanism. 



The second view is inspired by thermal infrared images of the atmospheres of Jupiter 
and Saturn at ~ 5 /mi (e.g, Westphal 1969; Westphal et al. 1974; Orton et al. 1996; Baines 
et al. 2005). In this spectral region gaseous opacity is very low, allowing flux from deeply 
seated, warm atmospheric regions to escape. Higher lying clouds, however, locally reduce the 
emergent flux. As a result these planets take on a mottled appearance with bright high-flux, 
low cloud opacity regions lying adjacent to cool er, darker, and cloudi er regions (optical depth 
varies from ~ to ~ 20 between these regions (jBanfield et al.lll998l )). lAckerman fc Marley 



(120011 ) suggested that the arrival of such cloud holes near the end of the L spectral sequence 
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may be responsible for the L to T transition. iBurgasser et al. I ( 120021 ) tested this hypothesis 
with a simple 'toy model' by summing weighted contributions of the spectra of cloudy and 
cloudless models. They showed that the observed J band brightening across the transition 
could arise from decreasing cloud coverage. Further support for this hypothesis comes from 
observations of L and T dwarf variab ility. While previous studies were somewhat equivocal 
(summarized in lArtigau et al.l (120091 )). two earl y T dwarfs have recently been shown to 
have large near-infrared photometric variability (lArtigau et al.l 120091 ; iRadigan et al.l |2010[ ) 
consistent with surface variations in cloud coverage modulated by rotation. 

The approach to modeling holes of Burgasser et al. was highly simplistic. The principal 
shortcoming being that it is not physically correct to combine the contributions of separate 
cloudy and cloudless models. Deep in the atmosphere of a brown dwarf the entropy in the 
convection zone must match that of the interior. Thus the temperature at a given, deep, 
pressure level is expected to be horizontally constant. However for a fixed T e g a cloudy 
atmosphere is everywhere hotter than a cloudless atmosphere. As an example Figure 1 
presents model a tmosphere profi l es for a un iform cloudy and a cloudles s atmosphere following 
the techniques of iMarley et al. I (120021 ) and ISaumon fc Marley I (120081 ). At depth, the cloudy 
profile is warmer than the cloudless profile by over 400 K; the difference in some models is 
even greater. Thus standard cloudy and cloudless models cannot simultaneously be valid 
descriptions of the real atmosphere at two locations even though both models are descriptions 
of an atmosphere with the same T e g. Clearly a new technique for self-consistently treating 
partly cloudy atmospheres is required. While we do not yet understand why cloud holes 
might appear, we here present a new approach inspired from models of Earth's atmosphere 
to model their influence and apply our results to model the spectra and colors of L and T 
dwarfs. 



2. Modeling Partly Cloudy Skies 

Instead of combining separately computed profiles for purely cloudy and cloud-free 
dwarfs we wish to construct a single, global temperature-pressure profile T(P) that incor- 
porates simultaneously the influences of both cloudy and cloud-free regions on the energy 
balance of the atmosphere. The final profile should conserve the total flux while allowing 
for nearby atmospheric regions to have differing cloud - but not thermal - profiles. This 
conceptually allows clouds to be displaced by winds, updrafts, or downdrafts and change 
location, as long as the global mean cloud fraction is constant. 

In three-dimensional terrestrial numerical weather prediction and climate models, clouds 
are typically smaller than the adopted computational grid scale. Various methods are used 
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to treat this in the radiative transfer calculations, one of which is to consider separate cloudy 
and cloud-free sub-columns in the same model column, with a single T(P) structure. This 
same approach can be used in a single column model of the Earth, allowing the global annual 
mean energy budget to be reproduced (Goldblatt & Zahnle 2010). 

Physically, such considerations are most important on Earth in the tropics. Taking the 
zonal mean T(P) and moisture profile would give a local runaway greenhouse due to the 
high water vapor content. However, dry areas in the tropics caused by subsiding air act as 
"radiator fins", allowing radiation from the surface to escape to space (Peirrehumbert 1995). 
In brown dwarfs atmospheres, cloud free areas would be analogous to the subsidence regions 
in Earth's tropics. 

To implement this approach we set the parameter < h < 1 to be the fraction of the 
atmosphere described by the cloud free sub-column, from which the local flux at some level in 
the atmosphere is J-hoic- The remaining fraction (1 — h) is the cloudy sub-column, with local 
flux Jcioud- Both sub-columns share the same T(P(z)) profile but have different emergent 
flux, in general J-hoie > -^cioud- We compute the total flux, J 7 t ot( z ) through the atmosphere, 
which is used in the radiative-convective equilibrium calculation, as 

Ttot{z) = hThoieiz) + (1 - /^cloudt^)- (1) 

We stress that this is not a combination of separate models, but rather conceptually repre- 
sents two adjacent sub-columns in the atmosphere with the same thermal profile and differing 
opacity. Our radiative-convective equilibrium model is then employed to solve for a single 
T(P) atmospheric thermal profile (McKay et al. 1989; Marley & McKay 1999) that carries 
net flux crT e 4 ff = J r to t(-2) through all statically stable layers of the atmosphere. Layers that 
would be convectively unstable in pure radiative-equilibrium are iteratively adjusted to fol- 
low an adiabat. The cloud profile and atmospheric chemistry are updated along with the 
atmospheric profile as it is converged to a solution, so the final thermal, chemical, and cloud 
profiles are all mutually self-consistent. 

A sampling of resulting model profiles is shown in Figure 1. The hottest profile is for 
a global, homogeneous / se( j = 2 model while the coolest model assumes no cloud opacity. 
Adding a small fraction of cloud free area to the f se d = 2 cloudy model efficiently cools the 
atmospheric profile. In the example shown setting h = 0.25 produces a model that lies over 
one-third of the way between the fully cloudy (/ se d = 2) and the cloudless (nc) models in the 
region above T{P) = T cff . Increasing the surface area fraction of holes to h = 0.5 produces 
an even cooler model. For comparison the figure also presents a homogeneous cloudy model 
with a much higher sedimentation efficiency, / se d = 4. This model with globally low cloud 
opacity follows a T(P) profile very similar to the partly cloudy model with thicker clouds 
(/sed = 2) with 50% of cloud-free regions. We have found that this is generally the case: a 
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partly cloudy model presents a similar thermal profile to a model with a thinner cloud, a 
point we return to in §4. 



3. Photometry and Spectra of Partly Cloudy Dwarfs 

The synthetic spectrum for a converged T(P) atmospheric profile is obtained from 



Equation ( 1). We compute abs o 



ute m agnitudes and colors with the ultracool dwarf evolution 



models of ISaumon Sz Marley I (120081 ). Figure 2 shows synthetic near-infrared photometry 
from the partly cloudly models. For a fixed T e Q, as h is increased the model colors briskly 
move to the blue in J — K and J — H . The J band flux increases as h increases from to 
0.75, but then dims slightly for cloud free models (h = 1). This is because atmospheres with 
even a small cloud cover are warmer in the atmospheric region from which the J band flux 
emerges. On the other hand, Mh is relatively constant across the transition from h = to 
h — 1 at constant T e g but shows a similar dimming as h — Y 1. 

The evolution of model dwarfs with fixed f se d produces trajectories that do not exhibit 
a rapid L to T transition as a global homogenous cloud sinks too grad ually below the photo- 



sphere (Burrows et al. 2006; Saumon & Marley 2008, Fig. 4). With the lAckerman &: Marley 



( 120011 ) cloud model, the L to T dwarf transition can only be modeled with an increase of 
the cloud sedimentation parameter / se< j. Figure 2 demonstrates that the transition can also 
be described as a progressive increase in cloud-free areas at a fixed / se d and T eS ~ 1200 K. 
The L/T transition dwarf colors and the J band brightening are well fit by this approach. 
However, 1200 K is slightly cooler than the observed T e g temperature of the transition of 
~ 1300 K (Golimowski et al. 2004; Stephens et al. 2009) and if T eS falls appreciably across 
the transition the J band bump would be weakened (Fig. 2). For a different choice of / se d 
and model gravity a different transition T c g would be expected. 

As noted in Saumon & Marley, there is an offset of ~ 0.3 to the blue in the J — K color 
of models from the bulk of L and T dwarf photometry. This may arise from shortcomings 
in the K band pressure-induced opacity of molecular hydrogen. Nonetheless, the trend of 
the late T dwarf J — K color can be better reproduced with models with h ~ 0.5 — 0.75. 
In the Mh vs. J — H diagram, the behavior of field L dwarfs and late T dwarfs are better 
reproduced by the cloudy / se d = 2 and the cloudless sequences, respectively although the 
latter would also be better matched with partly cloudy models with h ~ 0.5 — 0.75. This 
suggests the spectra of late T dwa rfs could be influenced by clouds, contrary to the usual 



assumption ( iBurgasser et al. Il2010f ). 
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4. Distinguishing Partly Cloudy Dwarfs 

In Saumon & Marley (2008) we demonstrated that by increasing the cloud sedimentation 
efficiency / se d as a brown dwarf cools from T e g = 1400 to 1200 K, the predicted model colors 
reproduce those across the L to T transition. In the previous section we likewise showed 
that increasing fractional cloudiness - at fixed / se d - has the same result. This leads us to 
consider how to distinguish the two cases. 

As shown in Figure 1 a partly cloudy T(P) profile (based on a / se d = 2 cloud) can be 
nearly identical to a model with a thinner homogeneous cloud but the spectrum from these 
two models are not necessarily the same because the former uses Eq. (1) to compute the 
flux. In the partly cloudy case some flux from deep, hot regions of the atmosphere (J\ \e 
is escaping through the clear regions that are otherwis e totally shielded by the clou d in the 



homogeneous case (see the middle panel of Figure 7 of lAckerman fc Marley I (120011 ) ) . Thus 



we expect that even for identical T(P) profiles the emission spectra will differ. 

Indeed that is the case as shown in Figure 4 which shows spectra computed from profiles 
shown in Figure 1. Focusing on the J band, which featu res the lowest molecular op acity 



and the deepest atmospheric window in the near-infrared ( lAckerman fc Marley 1 120011 ). the 
greatest flux is found for the cloudless model. The homogenous cloudy / se d = 2 model is 
faintest, with the / se d = 4 case falling in between. Even though the partly cloudy h = 0.5 
model has essentially the same thermal profile as the / se d = 4 model, it is brighter in 
the J band because some flux is escaping from deeper in the atmosphere. Of course flux 
conservation requires that the partly cloudy model must be fainter at other wavelengths, 
here in K band. Thus the model spectra are increasingly bluer in the near-infrared from 
/ sc d = 2 to / sed = 4 to the partly cloudy model. Figure 4 shows that the partly cloudy model 
is very close to the / se d = 4 homogeneous cloudy model at wavelengths where the flux is low 
and emitted from the upper atmosphere, and intermediate between / se d = 4 and cloudless 
in the JHK flux peaks where part of the flux comes from the deeper atmosphere. Thus, for 
a given T e g and gravity, the presence of holes in the cloud cover have a discernible effect on 
the near-infrared spectrum. 

For a given observed spectrum, where T e g, gravity and composition are not known a 
priori, how can we distinguish a partly from a homogeneous cloudy atmosphere? We explored 
this problem by fitting the near- infrar ed partly cloudy mod el spectra with our large library 



of cloudy models, using the method of ICushing et al.l (120081 ). and allowing T cff , gravity, and 



/ S cd to vary freely. We find that in general, the best fitting cloudy model has the same 
gravity, the same T e g (or slightly higher by ~ 100 K) and a higher / se d depending on h. The 
fitted cloudy spectra are close to the partly cloudy spectra, and some of the differences can 
be attributed to the grid spacing of the cloudy models. It appears that the JHK colors of 
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a partly cloudy model (T e g, g, / se( j, /i) can be well matched with a cloudy model (T e ' ff , 
/sed)- The differences in the near- infrared spectrum are thus subtle and com parable to the 



differe n ces found between obs erved spectra and the best fitting models (e.g. iCushing et al. 
(120081 ); IStephens et al. I (120091 ) ). Given that even with bolometric luminosity measurements 
T e ff temperatures cannot yet be measured to much better than ~ 50 K precision, effective 
temperature cannot yet provide a strong constraint on h. The most promising avenue to 
distinguish the two models in the L/T transition spectroscopically appears to be strong 
spectral features that are very sensitive to temperature and that are formed in the deeper, 
hotter regions of the atmospheres, such as the 0.99 jum band of FeH and the alkali doublets 
of Na I (1.14 /im) and K I (1.175 and 1.25 /im). Those should become stronger as the cloud 
coverage decreases. 

Given these considerations, variability may be the best method of characterizing cloud 
patchiness in brown dwarfs. Partly cloudy dwarfs may be variable if the lengthscales of the 
cloudy patches are large enough and if the viewing inclination is not near pole on. Variations 
in the geometry of the cloud cover ("weather") would also lead to detectable variations in 
brightness on time scales different from the rotation period. Indeed there is indication that 
a sizable fraction of L dwarfs are variable ( jGelino et al. 1120021 ; iKoen 1 120031 ) . Sample sizes 
are as yet too small to determine if blue L dwarfs or L/T transition objects show greater 
variability than other L dwarfs, but systematic studies would elucidate any trends. 

Polarization is another possible way to distinguish partly from fully cloudy objects. 
Clou d free brown dwarfs sh o uld n ot show appreciable polarization (jSengupta &: Marley 
20091 ). ISengupta fc Marley I (120101 ) demonstrated that very rapidly rotating, low gravity 
dwarfs with homogenous clouds can show sizeable polarization, but only for rotation periods 
less than about 4 hours. However partly cloudy L dw arfs may show a polarizati on signal at I 
band consistent with that observed in some objects (IMenard fc Delfosse 112004 ). Measurable 
polarization in a slowly rotating or high gravity brown dwarf would point to the presence of 
inhomogeneous cloud cover. 



5. Discussion and Conclusions 



There is growing evidence that the transiti on from the L to the T spectral class happens 
over a small effective temperature range (e.g., Stephens et al. Ifeood ). It is difficult for any 



model of a globally uniform, homogenous cloud to either sink out of sight or precipitate 
rapidly enough to account for the observed rapid change in spectral properties. Instead the 
transition may arise as holes appear in an otherwise uniform global cloud deck. 



-8 - 



We have presented the first self-consistent method for computing one- dimensional global 
atmospheric profiles appropriate to an atmosphere which has both clear and cloudy patches. 
The model does not depend upon the physical sizes of the patches, nor does it explain 
why patches might appear, but rather parameterizes the global fraction of clear and cloudy 
regions. The resulting T(P) profiles converged to by our model are intermediate in tem- 
pe rature between h o t, clou dy models and colder, cloudless ones. We confirm the conclusion 
of iBurgasser et al. I (120021 ). who used a more simplistic model to argue that patchiness in 
brown dwarf clouds is a plausible mechanism for the L to T type transition. 



In addition to the L to T transition, patchy clouds have been suggested by iFolkes et al. 



( 120071 ) as a po ssible mechanism to explain those L dwarfs that have unu s ually blue near- 



i nfrare d colors (ICruz et al.ll2003l ; iKnapp et al.l 12004 ; IBurgasser et al.l 120081 ). IBurgasser et al. 
( 120081 ) demonstrated that a T cff = 1700 K model with thin (/ sed = 4) clouds provides a 



good spectral match to the blue L dwarf 2MASS J1126-5003. We find that in this effective 
temperature regime a model with thicker ( f sef j = 2) and 50% cloud coverage produces a 



comparable spectrum. Thus we agree with IFolkes et al.l (120071 ) that partial cloudiness is 
potential alternative mechanism to thin global clouds for the blue L dwarfs. Differences 
in cloud coverage might result from different rotational rates, gravities, or even viewing 
geometries (pole vs. equator). 

The growing evidence for tempora l photometric variability at the L to T transition 
(lArtigau et al.ll2009l ; iRadigan et al.ll2010l ) supports the plausibility of partial cloudiness being 
responsible for the change in colors across the transitio n. Future 3-dimensional simulations 
of brown dwarf atmospheres (e.g., iFreytag et al.l |2010| ) coupled with a new generation of 
cloud models will likely help to elucidate the underlying mechanism responsible for cloud 
fragmentation at T e g- ~ 1200-1300 K. In the meantime additional observations of brown 
dwarf photometric variability and observations of L/T binaries will shed light on this poorly 
understood phase of their evolution. 
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Fig. 1. — Model atmosphere temperature-pressure profiles with T eff = 1400 K, logg(cm/s 2 ) = 
5 and solar metallicity. Solid lines are for atmospheres with horizontally homogenous cloud 
cover (labeled / se d = 2 and 4) or no clouds (labeled nc). Two partly cloudy models with 
hole fraction h = 0.25 and 0.5 and based on a / se( j = 2 cloud model are shown with dashed 
lines. The condensation curves of iron and forsterite are shown with dotted lines. 
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Fig. 2. — Near-infrared model color- magnitude diagrams (MKO system). Red and blue lines, 
respectively, show model colors for cloudy sequence with / se d = 2 and cloudless models (both 
logg = 5). Horizontal green lines connect partly cloudy model colors with T c g = 1000 to 
2000 K in steps of 200K (from bottom to top). Partly cloudy models (green dots) are for 
cloud-free fractions h of 0, 0.25, 0.5, 0.75 and 1, from right to left. Square symbols show field 
dwarfs with M dwarfs in bla c k, L d warf s in red and T dwar fs in blue. The photometry is 
pr imarily fromlLeggett et al.l (120021 ) and iKnapp et al.1 (120041 ) . Additional sources are given 
Saumon fc Marlev I (booi ) (Fig. 7). 
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Fig. 3. — Absolute model fluxes (d = 10 pc) for various cloud treatments, all with T e g = 
1400 K, log g = 5 and solar metallicity. Homogeneous cloudy spectra are shown for / se d = 2 
and 4 along with a cloudless spectrum. The partly cloudy model with h = 0.5 uses a / se( j = 2 
cloud. These spectra are computed from the thermal profiles shown in Figure 1. 



